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Abstract 

We propose a simple phenomenological model to estimate the spatial decoherence time in quan- 
. . tum dots. The dissipative phase space dynamics is described in terms of the density matrix and 

,J^ , the corresponding Wigner function, which are derived from a master equation with Lindblad oper- 
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ators linear in the canonical variables. The formalism was initially developed to describe diffusion 
and dissipation in deep inelastic heavy ion collisions, but also an application to quantum dots is 



i-^ . possible. It allows us to study the dependence of the decoherence rate on the dissipation strength, 

a ■ 

O \ the temperature and an external magnetic field, which is demonstrated in illustrative calculations 



on a circular GaAs one-electron quantum dot. 
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I. INTRODUCTION 

Decoherence processes in semiconductor quantum dots have attracted a lot of interest 
in the last years, not only due to their relevance for a quantum computer implementation^ 
but also because they present an experimentally accessible system to study the decoherence 
process in general^i"-. As demonstrated in several theoretical works^i""— , there are different 
processes that can lead to decoherence in a quantum dot, like interaction with optical and 
acoustic phonons or hyperfine interactions, in particular through electron spin coupling to 
a bath of nuclear spins. The processes occur on different time scales and are sensitive 
to external parameters like temperature or external magnetic fields. The coupling to the 
environment can be treated within the Born-Markov approximation^^ii^i^i, but also effects 
beyond the Markovian limit can play a role^^'—'^'^. However, the vast majority of the 
processes studied so far focuses on spin decoherence, mainly because it is the spin of the 
electron(s) that makes a quantum computational application of quantum dots possible^. 
Nevertheless, also spatial decoherence which is arising from dissipative phase space dynamics 
in the canonically conjugated coordinates and momenta (see e.g.— for a detailed discussion), 
can possibly become relevant. As an example, one can think of a situation where the space- 
and spin part of a wave function are connected by the fermionic total asymmetry condition. 
Also in such a scheme as recently suggested in^-, where electromagnetic transitions in solid 
state devices are used for controlled operations, phase space dynamics can be important 
since it is not only the spin but also the total angular momentum that plays a crucial role. 

In the present work, we aim to estimate the spatial decoherence time scale and to study 
its dependence on the coupling strength to the environment, on the temperature and on an 
external magnetic field. The latter plays a role since it determines the cyclotron frequency 
and hence also the effective confinement strength, which, together with the temperature, 
was shown to influence the asymptotic spatial decoherence in quadratic potentials^. In 
addition, the magnetic field also explicitly influences the time evolution of the system. 

The study is carried out using an analytical model in the Markovian limit with linear 
Lindblad operators, which was initially developed to study diffusion and dissipation in heavy 
ion collisions^i^i. We show that, with an appropriate choice of the involved constants, the 
model can also be used to describe a two-dimensional quantum dot in a perpendicularly 
applied external magnetic field (section [TTl)- The temperature dependence is incorporated 



in the diffusion coefficients which significantly determine the time behavior of the density 
matrix^. The model has the advantage that it is rather general and allows us to include 
environmental effects without the necessity to explicitly compute the system-environment- 
interaction. The latter is instead taken into account by phenomenological constants that 
emerge from the Lindblad operators. This is sometimes referred to as the reduced dynamics 
approach. In practical calculations, it is then only necessary to find appropriate values for 
these constants, depending on the environmental effects one wishes to consider. In the illus- 
trative calculations shown here, we relate these effects to the electron-phonon-interaction, 
but one could also consider other effects within basically the same model without loss of 
generality, provided the required input (as, e.g., electron-phonon scattering rates in our case) 
is known at least approximately. In section IIIIl we discuss how the decoherence time scale 
can be determined from the model, using the general results derived in Ref.— . It is demon- 
strated in more detail in illustrative calculations on a circular GaAs one-electron quantum 
dot in section IIVI followed by a discussion and concluding remarks in section |Vl 

II. DESCRIPTION OF THE MODEL 

The Hamiltonian of a one-electron quantum dot with a harmonic confinement (which is 
a very common choice, see e.g.—"—) in the xy-plane and exposed to an external magnetic 
field B pointing in z-direction reads 

2m* 2 ° 8m* 2m* ^ ^ ^ ^ 

Here, uq denotes the confinement strength, r = a/x^ + y"^ is the radial polar coordinate, L^ 
and Sz denote the 2;-component of the orbital angular momentum operator and spin operator, 
respectively, fit, is the Bohr magneton and m*, g* are the effective mass and effective (/-factor 
for the used semiconductor. By defining the cyclotron frequency Uc = eB/m* and the 
effective frequency u = \^uiq + lu^/A, it can be written as 

H = Ho + Hs (2) 

where Hg = g*fibBSz is the spin part and 
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The latter expression is a general Hamiltonian of a two-dimensional harmonic oscillator 
exposed to a perpendicular magnetic field. The phase space dynamics of such a system 
when coupled to the environment was, for example, studied semiclassically in^ by means of 
a Fokker-Planck equation. Also an explicit inclusion of a heat bath in this Hamiltonian is 
possible, which leads to non-Markovian dynamics as recently demonstrated in^I for a nuclear 
system. Another possible approach is the infiuence functional method'^^, which was used to 
study the decoherence dynamics of two coupled harmonic oscillators in a general environ- 
ment, with their potential minima being separated by a finite distance^^. Here, however, we 
will adopt a simpler phenomenological picture to describe the quantum dot coupling to an 
external environment, based on a Markovian master equation^ which is known to be valid 
in the weak coupling limit^. In this framework, dissipation and decoherence are described 
by Lindblad operators, which is a rather common approacb^Sii^^ii. The formalism is based 
on earlier work of Gupta et a^ and Sandulescu et a?^ which is briefiy introduced below. 

Due to excitation of internal degrees of freedom (i.e. the nucleons) in heavy ion collisions, 
dissipation is a rather important issue in its quantum mechanical description. It is common 
to describe such a dynamics in terms of dimensionless coordinates of proton and neutron 
asymmetry, defined as gi = {Zi — Z2)/{Zi + Z2), q2 = (^1 — ^2)/(^i + ^2) where Zi, Z2 
and Ai,A2 are the charges and masses of the colliding nuclei. A model to couple these 
coordinates was suggested in^-. Later, this model was generalized io^, where the complete 
description of the dissipative dynamics was explicitly derived from the Markovian master 
equation for the density matrix p given below, 

j 

where Vj is a set of Lindblad operators, and the considered Hamiltonian had the following 
form: 

El L 2 ^^k^k 2 1 \ — ^ 

( ^^Pk + -^—1k ) + 2 2^ ^^I'^a (Pfcife + fePfci) + i'i2qiq2 + K12P1P2. (5) 
fc=l ^ '^ ^ fci,A:2 = l 

Here, pi,P2 are the canonically conjugated momenta to the charge and mass asymmetry 
coordinates. The appearing coupling constants can be partly calculated from the nuclear 
liquid drop model or determined by fitting to experimental data. However, if these constants 
are chosen as z/12 = 0, kx2 = 0, /in = 0, /i22 = and 

qi = X, q2 = y, pi = Px, P2 = Py, mi = m = m2, Ui = u = uj2, /i2i = ^-^ = -pu, (6) 



the Hamiltonian i^o from Eq. ([3]) is reproduced exactly. Moreover, following the same choice 
of linear Lindblad operators as in^i, 
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Vj = J^ (ajkPk + bjkQk) , j = 1, 2, 3, 4 (7) 

fe=i 

where ajk, bjk are complex numbers, we see that the spin part Hg of the full quantum dot 
Hamiltonian ([2]) commutes with the full Hamiltonian as well as with the Lindblad operators, 
so that the resulting equations of motion of the first and second moments in the canonical 
variables are unaffected. In other words, all results derived in^ also remain valid in our 
case. Thus, we will omit the derivation, solution and discussion of the equations of motion 
here and only briefly quote the main results relevant for our study. At this point, we would 
also like to mention that, since the spin motion decouples, no spin dephasing effects are 
present in the spatial decoherence studied here. The spin relaxation and dephasing times 
(often called Ti and T2) can also be studied within an equations-of-motion approach (see 
e.g. Ref.— and the references within), which, however, requires different models and is not 
considered in the present work. 
We use the abbreviations 



mit) = {{x),{y),{p^),{py)f 
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for the time-dependent expectation values of the canonical phase space operators and 
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for the time-dependent symmetric covariance matrix, where the elements are defined as 



(TAB = OBA = ^(^^ + ^^) - (^)(^) 



(10) 



for any two operators A, B. The time evolution of the expectation values is given by 



m(t) = exp(ty)m(0). 



(11) 



where m(0) denotes the initial expectation values and Y is the time evolution matrix, which. 



after insertion of the constants given in ([6]) into the general result from^, becomes 
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The phenomenological dissipation constants Xm, ^12 and /3i2 emerge from the Lindblad op- 
erators ([7]) and are exphcitly given by 



ai2 = -Im(ai,a2), 
/3i2 = -Im(bi,b2), 
Xki = -Im(ak,bi), 

where the vectors ak,bi are defined as (cf. Eq. ([7])) 



(13) 



^k — (ctlfc, 0'2k, Osfc, 0"ik) , bk — {bik, &2A;, ^Sfc, &4fc) 



(14) 



with the scalar product 



(x,y) = E 
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However, since we consider a circular one electron quantum dot where the dynamics is 
symmetric in x and y, we set the off-diagonal dissipation constants to be zero here and in the 
following (i.e. A12 = = A21 and ai2 = = (3i2) and, furthermore, demand An = A22 = A, 
hereby restricting the dissipation strength to a single phenomenological parameter. For the 
covariance matrix, the following time evolution is derived: 



a{t) = exp(tF)(a(0) - a{oo)){exp{tY))'^ + a{oo). 



(16) 



Here, cr(0) is the initial covariance matrix and a (00) its asymptote. The latter can be 
determined from a set of diffusion coefficients, which are given by 
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where the notation gi = x,q2 = y,pi = Px,P2 = Vy is used. They are connected to the 
asymptotic covariance matrix by the relation 



Yaioo) + a{oo)Y' = -2D, 



(18) 



where D is the symmetric diffusion matrix 
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The choice of the diffusion coefficients is, in general, a non-trivial issue, since there are 
several conditions that have to be obeyed in order to preserve the non-negativity of the 
density matrix and the uncertainity relation. This will not be discussed here (see Q.<gM-^— 
for more details). In the present work, we use a two-dimensional extension of the commonly 
used temperature-dependent coefficients of a harmonic oscillator without further mixing, 
such that the diffusion matrix is diagonal: 
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and Dij = otherwise, where T is the temperature and k the Boltzmann constant. From the 
given time evolution of the first and second moments, one can obtain the Wigner function fw 
of the system, which is the best possible quantum mechanical analogon to a classical phase 
space density (although it is, in general, not positive everywhere and, therefore, cannot be 
interpreted as a true density). The latter was found by means of Weyl operators in^: 

1 f 1 
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where ^ = {x,y,px,Py)^ is the phase space vector. This agrees with the result obtained in 
earlier work— i^*^. In the following, we use this result to calculate the decoherence rate. 



III. DECOHERENCE RATE 



It is by far not trivial to give a general definition of 'decoherence'. Often, it is simply 
referred to as 'loss of coherence' in a quantum system or as entanglement of the latter with 




its envitonment. Technically, however, the degree of decoherence can be expressed through 
the density matrix of a quantum system, or, more precisely, through the damping of its 
off-diagonal elements^'^'^"— . We shall adopt this definition in the following, although it is 
also possible to study decoherence directly in terms of the Wigner function^^. The density 
matrix is connected to the Wigner function by the transformation 

{x,y\p\x\y'){t)= (22) 

dp^dpy exp {-{px{x - x') + py{y - y'))\ fwiix + x')/2,{y + y')/2,p^,py,t), 

which can be carried out analytically for the Wigner function (!2T|) . After evaluating the 
two-dimensional Gaussian integral and introducing new coordinates 

_ x + x' . _ _ , 

S. = ^, Ay = y-y\ (23) 

we arrive at the following expression: 

{x,y\p\x',y'){t) = AT exp [-K^iE, - {x)Y - K^iJ^y - {y)f - K,{^x - (x))(S, - {y))- 
K^Al - K,Al + K^KAy + iK-j{^^ - (a;))A, + iK^{^y - (y)) A,+ 

(24) 
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The normalization factor is given by 
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where Cij = Cji denote the elements of the inverse of the covariance matrix cr(t) ^. The 
explicit form of the time-dependent constants Ki{t)..Kio{t) in terms of the matrix elements 



Cij{t) is given below: 
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Note that also the expectation values are functions of time, governed by ( TTTl) . In the par- 
ticular case of a circular one electron quantum dot, the dynamics is considerably simplified 
due to the symmetry in x and y, because the dispersions of the terms quadratic in E^., T,y 
that determine the damping of the diagonal elements of the density matrix are identical 
{Ki = K2). The same holds for the dispersions of the terms quadratic in A^^, Aj^ that de- 
scribe the damping of the off-diagonal elements (7^4 = -ft's), which, foUowing^^i^"— , allows 
us to define a single decoherence parameter 

The definition is such that (5/5 = 1 corresponds to a perfectly coherent state and 60 = 
implies that coherence is lost. In the next section, we investigate its time behaviour under 
the influence of the magnetic field, the temperature and dissipation. At this point, we 
would like to emphasize that not only the decoherence degree but also other purely quantum 
mechanical quantities can be extracted from the Wigner function in the case studied here. If 
the quantum dot is prepared in a coherent state, the Wigner function is positive everywhere, 
and hence it coincides with the classical phase space density^. For example, the quantum 



mechanical probability density p(x, y, t) is directly obtained from the Wigner function by 
setting X = x', y = y' (i.e. T.^ = x, T,y = y, A^ = = Ay) in Eq. fl24|) . 

IV. ILLUSTRATIVE CALCULATIONS 

As an illustration, we consider a one-electron GaAs-quantum dot. Throughout this sec- 
tion, we use effective atomic units, i.e. atomic units scaled with the GaAs material param- 
eters er = 12 A (dielectric constant) and m* = 0.067 Trig where rrie is the electron mass, and 
a confinement strength of hwQ = 5meV. Disregarding the spin part, which, as previously 
mentioned, is irrelevant for the phase space dynamics studied here, the solutions to the 
Hamiltonian ([3]) can be written in terms of a principal and an angular quantum number 

^nm{r,^) = Unm{r)e'"'''- (28) 

Since we consider the system to be initially in a prepared state, the initial conditions required 
for the time evolution of the first and second moments (Eqs. [11] and [T6l) are, unless given 
analytically, calculated numerically for any state (1251) using a B-Spline basis^ as demon- 
strated in-. In particular, we have m(0) = for the expectation values of x,y,Px,Py and 
the initial covariance matrix has the following form: 
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where a^x = <^yy and Op^p^ = Cpypy (the equality is again a consequence of the symmetry) 
are numerically calculated values. They depend on the effective confinement frequency 
and therefore on the magnetic field since it influences the latter {u = a/wq + i^c/4 where 
cOc = eB/m,*, cf. section [TTl)- The explicit dependence is illustrated in figure [H This 
behavior can be understood qualitatively by considering a usual one-dimensional harmonic 
oscillator, where the variances of the n-th state are given by a^x = (1 + 2n)h/{2m*uj) and 
'^PxPx = (1 + 2n)hm*uj/2. Physically, this simply means that with increasing magnetic 
field the electron becomes more localized in position space and less localized in momentum 
space. The uncertainty relation in each coordinate, however, does not depend on uj (and 
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FIG. 1: (Color online) Variances in position and momentum shown for some one-electron GaAs 
quantum dot states \nm), given in effective atomic units: Uxx has dimension length^, where length 
is measured in effective Bohr radii (oq ~ 9.794 nm) and CTp^p^ is measured in Hartree*m,* where 
m* = 0.067 me and Hartree* w 11.857 meV is the effective atomic energy unit. 

hence neither on the magnetic field) since it is determined by the product of the variances: 

1 






(30) 



(and in the same way for y,Py). Also, for the quantum dot states, the initial covarinaces 
a^p^ and ayp^ always vanish. Among all states, it is only the ground state \nm) = |0 0) 
that has minimum uncertainty due to its ideal Gaussian shape, while for excited states the 
Gaussian shape is disturbed and the uncertainty relation becomes a strict inequality. This 
corresponds to the fact that the ground state of a harmonic oscillator is a particular case 
of a Glauber coherent state. Hence, ^^(t = 0) must be equal to unity for the quantum dot 
ground state which agrees with the calculations and retrospectively confirms the imposed 
definition of the decoherence parameter. 

In addition to the previously discussed initial conditions, also the coupling strength to 
the environment needs to be determined. As already mentioned in the introduction, the ad- 
vantage of the model used here is that one can phenomenologically account for, in principle, 
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FIG. 2: Time evolution of the probability density in coordinate space shown for B = 1T,T = 25K 
and A = lO^'^w where the dot is initially prepared in the ground state. The coordinates are given 
in effective Bohr radii Gq ~ 9.794 nm and the probability density in units of l/(aQ)'^. 

any kind of environmental effects simply by choosing appropriate values of the phenomeno- 
logical constants. However, since the master equation (jl]) is only valid for weak coupling of 
the reduced system to the environment, the ratio X/u should be much smaller than unity, 
which restricts the model to the description of processes that obey this condition. Here, we 
consider the dissipative effects to be caused by electron-phonon interactions and therefore 
the dissipation rate A is approximately given by the electron-phonon scattering rate. The 
latter ones were found to be typically of the order A ~ 10 ns~^ in GaAs structures^"—. For 
the effective confinment strength chosen here, this yields A ~ lO^'^cj and hence the Marko- 
vian approach can be seen as justified in the present study. Another aspect to be considered 
is that, apart from being incorporated in the diffusion coefficients (Eq. |20|) . the temperature 
also influences the phonon scattering rate. However, since the phonon scattering rates for 
different temperatures are not known exactly, the parameter A is varied independently in 
the relevant region at different temperatures to investigate its influence on the decoherence 
time scale. 
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FIG. 3: (Color online) Time evolution of the decoherence degree for three different environmental 
coupling strengths A. The solid curves correspond to a temperature of T = 25 K, the dashed curves 
to T = 50K and the dotted curves to T = 100 K, while the different magnetic fields are indicated 
by symbols: i? = 1 T (circles), B = 3T (squares) and B = 7T (triangles). 

Figure |2] shows the time evolution of the probability density (obtained by setting x = 
x', y = y' in Eq. f l24|) ) and figure [3] the time evolution of the decoherence degree when the 
quantum dot is initially prepared in the ground state. We observe that the asymptotic degree 
of decoherence strongly depends on the temperature and, more weakly, on the magnetic 
field strength. Naturally, higher temperatures lead to stronger decoherence, while a strong 
magnetic field has a protective effect on quantum coherence. This corresponds to the result 
obtained for the one-dimensional case^^, where the asymptotic decoherence degree was shown 
to be 6£,{oo) = ta.nh.{hw / {2kT)) (Recall that the effective confinement u increases with B). 
However, the time scale itself on which the system approaches its asymptotical decoherence 
is determined by the coupling to the environment. Depending on the latter, decoherence 
occurs on a time scale between 20 — 200 ps. It is interesting to note that the rule of thumb 
to estimate the ratio of relaxation (r^) and decoherence (td) time scaled 
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where LdB is the thermal de Broghe wave length, 



MB 



(32) 



\/2m*kT 

holds in the case studied here, but the ratio is of the order of unity, which is a rather untypical 
behavior. In fact, in many situations decoherence is several orders of magnitude faster than 
dissipation and relaxation; For macroscopic systems, the ratio ( 13T|) becomes astronomically 
large^>^>^. In quantum dots, however, the wave packet spread can be of the same order 
as the thermal de Broglie wave length - a very remarkable property, once again displaying 
their fascinating features. 

V. CONCLUSIONS 

Using a Markovian master equation approach with linear Lindblad operators, we inves- 
tigated the dissipative phase space dynamics of a one-electron quantum dot. We obtained 
the density matrix in coordinate representation from which an expression for the spatial 
decoherence parameter was derived. With numerically calculated initial values for the first 
and second moments of a quantum dot in a prepared state, we analyzed the time evolution 
of decoherence and the influence of the temperature and an external magnetic field. The 
phenomenological coupling to the environment was assumed to emerge from electron-phonon 
scattering. We found that the asymptotic decoherence strongly depends on the temperature 
and also on the magnetic field and the decoherence time scale is driven by environmental 
coupling strength, varying between a few and a few hundred picoseconds. 

The model presented here has the advantage of being rather simple. The price one 
has to pay is, on the one hand, its phenomenological nature and, on the other hand, the 
restriction to Markovian dynamics. Hence, the obtained results should be viewed with these 
limitations in mind. If, for example, one aims to investigate decoherence arising from faster 
processes than electron-phonon scattering, the weak coupling limit is not necessarily valid, 
since the dissipation constant approaches the confinement strength. At the same time, the 
phenomenological nature also has the advantage that one is able to describe different kinds 
of interactions without loss of generality, just by choosing different values for the emerging 
constants, as long as the Markovian condition is not violated. 

Further difficulties can occur, if, for example, the circular symmetry is disturbed or an- 
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harmonic effects have to be accounted for. In that case, the off-diagonal couphng constants 
in Eqs. flT2|) or flT9|) may be non-zero, and it may also be quite non-trivial to find a sin- 
gle decoherence parameter in this case. It should be stressed that analytical solvability is 
provided for a purely harmonic confinement only, while for more complex potentials even a 
numerical solution is not always possible, since no finite system of the equations of motions 
can be derived. 

Another conclusion we can draw from the present study is that, if the temperature is very 
much smaller than the confinement strength {kT ^ hu), spatial decoherence should not be of 
significant relevance. Even at a temperature of 25 K we see that the asymptotic decoherence 
is about 6r)(t — )■ oo) ~ 0.9. However, in some experimental setups the temperatures are as 
low as T = 100 mK, and hence we can conclude that spatial decoherence practically does 
not occur in that case. 

As for further applications of the model presented here, it should be also suited to study 
tunneling processes in quantum dot molecules. Calculations on tunneling through one- 
dimensional parabolic potentials with dissipation described by linear Lindblad operators 
were, e.g., demonstrated in^*^, which can be extended to the two-dimensional case. 
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